Introduction to Quantum Field Theory 


Daniel N. Blaschke 


Faculty of Physics, University of Vienna, 
Boltzmanngasse 5 A-1090 Vienna (Austria) 


Version: May 2, 2011 


Foreword 


These lecture notes are meant to give a concise introduction to Quantum Field Theory 
using path integral methods. They are mainly based on the books by W. Cottingham and 
D. Greenwood [I] (Chapter 1), A. Das [2] (Chapter 2) and L. Ryder [3] (Chapter 3) with 
supplemental material from Refs. [41/6]. The interested reader may want to consult these 
books for further details. 


© D.N. Blaschke, May 2, 2011 
E-mail: daniel.blaschke@univie.ac.at 


Contents 


CL. Introduction 4 
1.2.1. From particle mechanics to field theory ea oe aoe ee eee S) 
7 
8 
9 
0 


[1.2.3. Scalar field (heard... 
(1.2.4. Dirac fields 


11.2. The role of symmetries in physicd Ue ce de ee ge a 1 
2. Path integral methods 12 


2.1. Path integral formulation in quantum mechanicg . . . oa o a aa aaa 12 


[2.2. Generating Functional] Beene op BTN ae ee ee tee Ee, tte et Bg aad a 13 
2.3. Path integral for Grassmann variables ..........0.0.00 000020 a se 14 


[2.4. Path integral for a relativistic scalar field theory| ne oS ae 8 ete, HS ha ae 15 
[2.5 Effective action! 


[3.1. Regularization and power counting eek te oe oe he eee ee eee a ee ee 29 
3.2. Renormalization and renormalization group... 2... 6. ee ee 31 


A. Supplemental Material 34 
[A.L. Dirac formalism 


1. Introduction 


1.1. A brief history of Quantum Field Theory 


Quantum electrodynamics was formulated about 1950, many years after Planck’s original 
hypothesis (1901) that the electromagnetic field should be quantized. In subsequent years, 
quantum field theories were studied extensively and one was finally lead to the “Standard 
Model” of particle physics which describes elementary particles and their interactions in flat 
space-time. Three such basic interactions are involved: 


1. the electromagnetic field, whose gauge bosons are massless spin 1 “photons”, 


2. the weak field, whose gauge bosons are the massive spin 1 WF, Z bosons, 


3. and the strong field (QCD), whose gauge bosons are the (postulated) 8 massless spin 
1 “gluons”. 


The first two forces are described by the Weinberg-Salam electroweak model based on the 
gauge group U(1) x SU(2), and the W*, Z bosons become massive due to spontaneous 
symmetry breaking through the so-called Higgs mechanism. The strong force, on the other 
hand is based on the gauge group SU (3) whose three (anti)charges are dubbed (anti)red, 
(anti)green and (anti)blue, hence the term “chromo” (colour) in the models name. 

Additionally, the elementary particles charged under these forces can be categorized into 
three generations of quark-doublets and equally many generations of lepton doublets. Details 
are given in the Table 


Table 1.1.: The elementary particle zoo 


1. generation 2. gen. 3. gen. 


a| p [down (a) [| charm (o) [strange ©) top (1) [bottom O) 
charge 2 
[nas er a REEE SO WOME] 


lepton] electron e~ | el.-neutrino Ve] muon p | p-neutr. U | ae T-neutr. Vr 


charge] -1 | 0 


E ain | R | toner arar [= sone] 


Many experiments, mostly particle accelerators, have verified the validity of the standard 
model. The only unobserved particle to date is the Higgs particle for which the search is 
ongoing. Furthermore, the standard model has known shortcomings: In its original form, 
massless neutrinos are considered which meanwhile is known to be untrue. Neutrinos do 
in fact have mass and according extensions to the standard model have been suggested. 
Furthermore, one expects that at energies at/or beyond the TeV scale new physical effects 
which are not covered by the standard model will appear. 

Some important discoveries were: 


e the electron e~ 1897 by J.J. Thomson, 
e the anti-electron (or positron e*) 1932 by C.D. Anderson in cosmic radiation, 
e and in the same year J. Chadwick discovered the neutron, 


s the muon ~~ 1936 by C.D. Anderson in cosmic rays, 


e the WF, Z bosons (80/90GeV) in 1983 at CERN, 


just to name a few. 


History and overview of some particle accelerators. 


1. Early 1930s: Cockroft-Walton linear accelerator at Cambridge, UK (0.7 MeV); and 
Lawrence’s cyclotron at Berkley, USA (1.2 MeV); 


2. further accelerators in the 1950s and after; 

3. TEVATRON (Fermilab), p-p collisions, 900 GeV, 1987-2011, discovered top-quark; 
4. SLC (SLAC, Stanford), e~-e* collisions, 50 GeV, 1989-1998; 

5. HERA (DESY, Hamburg) e at 30 GeV, p at 820 GeV, 1992-2007; 

6. LEP2 (CERN), e~, et, 81 GeV, 1996-2000; 


7. Vienna Environmental Research Accelerator (VERA), 3 MV tandem accelerator, for 
Accelerator Mass Spectrometry (AMS), 1996-present; 


8. PEP-II (SLAC, Stanford), e~, e+, 9 GeV, 1999-2008, 


9. Relativistic Heavy Ion Collider “RHIC” (Brookhaven National Laboratory, New York), 
2000-present; 


10. LHC (CERN), p, 7 TeV, 2009-present. 


Concepts of renormalization/regularization. 

Charged particles receive a “self-interaction”, i.e. they feel their own field. “Naked 
charges” can never be measured, only the effective charge including self-interaction of the 
particle. In fact, the value is energy dependent (RG-flow) as described by the so-called 
6-function. For example, at low energies the fine-structure constant of QED is given by 

e? 1 


= x — 1.1 
9 Areghc 187’ cit) 


but at high energies a becomes larger, as can be measured in particle accelerators. This can 
be interpreted by the following gedanken-experiment: vacuum fluctuations constantly lead 
to the production and annihilation of (virtual) particles and antiparticles. Near a charged 
particle, say an electron, these will shield some portion of the charge leading e.g. to the 
value of 1/137 at low energies. Measuring at high energies means we can detect the charge of 
that electron at closer distance and less shielding, hence the increased value. In fact, being 


interpreted as a point particle, the beta-function diverges as the energy goes to infinity. 
However, this picture is not expected to hold at arbitrarily high energies. 

On the other hand, there are also quantum field theories, where the beta-function behaves 
in the opposite way: In QCD, for example, the coupling is strong at low energies and becomes 
weak at high energies. (Such theories are called asymptotically free.) The mathematical 
reason is that quantumchromodynamics QCD is a non-Abelian gauge theory, meaning its 
gauge bosons are charged and interact with each other. Therefore, the naive “screening” 
picture above does not apply here. In fact, so-called “confinement” prevents production of 
charged particles. 

These ideas will be made more explicit in Section 


Notation. 
Throughout these lecture notes, natural units as described in Appendix|A.2] are used. 


1.2. Classical theory 


1.2.1. From particle mechanics to field theory 


Consider a point particle of mass m whose position at time t is x(t), and which is driven by 
a force F = V' (x). Its equation of motion (i.e. Newton’s second law) is then given by 


dx dV 
— L n =0. 1.2 
MH t dz ý (1.2) 


A way of deriving this equation is by the principle of least action. The action S' in the above 
case is given by 


t2 
m (dx\? 
s= lar E SPO (1.3) 
tı 


where L is called the Lagrangian, and the integral in the action is taken over a path from 
tı to t2. The principal of least action states that the path the particle actually takes (from 
the infinitely many possible ones) is the one for which S is a minimum. It can be derived 
by varying S and requiring this variation to vanish, i.e. 


t2 

dx dôx dV 

œl (EE Tn (1.4) 
ty 


Integrating the first term by parts keeping in mind that the variation is zero at the end 


points t1,2, we find 
fi dr dV 
Sa Ea = 1. 
lal Te 0, (1.5) 


which for arbitrary variations 6x implies Eqn. (1.2) above. The above can of course readily be 
extended to a system of particles, or a more complicated mechanical problem with n degrees 


of freedom q;. The Lagrangian then becomes L(q;, ġi). If we push this generalization further, 
say to an infinite number of degrees of freedom, we replace the q(t) with fields ¢(x,) = 
o(t,x,y,z) where the index i is replaced by continuous variables x,y,z. Furthermore, we 
define the Lagrangian density £ by L = f@al, and a general field theory action hence is 
written as 


g= J PENO) (1.6) 


depending on fields various fields ¢; (which might have additional indices) and their gradi- 
ents. We furthermore assume “natural boundary conditions”, i.e. that these fields vanish 
sufficiently fast at infinity to justify integration by parts. Variation of this action leads to 


OL OL 
as = fate (SS - aoe) 66%, 1.7 
op DO) °° aA 
implying the generalized Euler-Lagrange equations 
L ee =0. (1.8) 


06; "O(n di) 


1.2.2. Maxwell fields 


In the case of photons, the action is given by 
1 
S= [es (irer +u") , (1.9) 


where A,, is the vector potential, j, denotes an (external) current density and 


ù -E =E =E 
E, 0 -B, B; 
Ey Be 0 -B |’ 
E, —By Bz 0 


F” = OH AY — AV A! = (1.10) 


is the electromagnetic field strength tensor in vacuum. Furthermore, we use the Lorentz 
metric Nav = (+, —, —, —) for pulling indices up and down, and consider natural units where 
es h= Lick Appendix[A.2). Varying the action with respect to the A,, leads to the 
equations of motion] 


ð, FP” = 0,0" AF — o” (ð, A”) = =a Ñ (1.11) 


These are in fact the inhomogeneous Maxwell equations, i.e. for u = 0 one has with 


> 


B=-VAo—2 and p = j°: 
divE = p, (1.12) 
whereas using B =rotA the U = 1, 2,3 components yield 
>- OR >- 
LB =j. 1.13 
ro a (1.13) 


‘Note, that these imply current conservation since 0 = 0,0, FH” = ð j". 


The homogeneous Maxwell equations follow as an identity directly from the definition of 
the antisymmetric field strength tensor (1.10), 1.e. 


~ 1 
oF” = Ou Fon See" 0.0. 4. = 0. (1.14) 


The u = 0 and u = 1,2,3 components correspond to 


divB=0 and roat E + = =0, (1.15) 
respectively. 


Exercise 1 Derive the Maxwell equations starting from (1.11) and (1.14). 


The action (1.9) as well as the field strength tensor F» are invariant under the “gauge 
transformation” 


Al, = Au + Oar, (1.16) 


(up to a surface term provided the current density is conserved, ô j” = 0). This gauge 
freedom can be fixed by demanding a gauge condition, such as e.g. the Lorenz condition 
ð, A" = 0. 


1.2.3. Scalar field theory 


Let us now consider the simplest field theory, which represents a good toy model in order to 
learn general properties of quantum field theory. We start with a free field whose classical 
action is given by 


S4 = [as (-3 „pao + srg? = io) (1.17) 


The e.o.m. follow as 


(+m?) (2) = j(2). (1.18) 


A solution can then be found by employing the Green function method where 


glz) = 0%) + lage, dic). 
(0 +m?) Glz, x) = 64(z —2'), (1.19) 


where $° solves the homogeneous equation and is chosen in such a way that 0 satisfies the 
boundary conditions. Making use of translational invariance, the first line of (1.19) can be 
replaced by an algebraic one by making a Fourier transformation leading to 


dtp cp) 


G(x —2') = — oe (1.20) 


Crip mè 


Obviously one has to deal with the zero of the expression pê — p? — m? by prescribing a 
slightly deformed contour of integration. This in turn leads to the definition of the retarded 
and advanced Green functions, as follows: 


d£ ip(z—z2’) 1 1 
G ret (a T x’) = -J p is -M 
adv (27)* 2wp Po—wptie po+wpie 
( 


a — y! 3 ser t 00 
_ BL = 2) [Pere cin (wp(zo — 24) , (1.21) 


— 


T 
£ 
2T)? Wp 
where wp = y P? + m?. 


In order to find the homogeneous solution ¢° we once more employ a Fourier ansatz 
leading to 


4 
oa | 55 SALE — m?)¢"(p) 
S 3 . 
~ a T we (ot (P) + ope”) , (1.22) 


where ¢*(p) := $'(+wp, £P). 


1.2.4. Dirac fields 


The Schrödinger equation for e.g. an electron wave function 10 is given by 


OW 
i1— = 


a TEY. (1.23) 


In order to ensure a symmetry between space and time (as required by special relativity), 
Dirac postulated the Hamiltonian for a free electron to be of the form 


H = -igV + Bm, (1.24) 
so that 
ð e 
iz + iav — pm v= 0, (1.25) 


where a; and 8 must be matrices satisfying a Clifford algebra, since the solutions to the 
above equation should also be solutions to the Klein-Gordon equation leading to 


L =iay + Bm) (iz + IY — Bm) w 


07 1 ; 
= -f +3 {ai aj} 0,0; + im {04, B} 0; — Bm?) Y 
=0, (1.26) 
where 
{ai aj} = 2ðij , B? — 1, {a;, B} =0. (1.27) 


Though not unique, a possible representation is the so-called chiral representation where 


_ [(-c 0 Mhk 
e TE C na 
and ø are the usual Pauli matrices. We may now rewrite Eqn. (1.25) as 
o > 
(ið L- m) p= B (5 + iaV — Bm) p=0, (1.29) 


where y° = 8 and 7f = Ba; and {y",y”} = 2n”, u,v € {0,1,2,3}. The Lagrangian must 
be Hermitian and scalar leading to the above e.o.m., and hence 


£ Bi", - m) Y, (1.30) 


where Y% = w'y7°. Another useful matrix is y° = iy°y!y?73 which in the chiral representation 


is given by 
-l2 U 
5 2 
y =( 0 4 S (1.31) 


The matrices 4 (14 Æ 7) are projectors giving the right and left handed parts of a Dirac 
spinor, i.e. 


nmu UE).  ZQatr)v=()). (1.32) 


The Dirac equation describes particles with intrinsic angular momentum Žo and intrinsic 
magnetic moment Lo if the particle carries charge q. Furthermore, there exist negative 
energy solutions which Dirac interpreted as antiparticles. This, of course at the time was 
a rather bold claim, which however was experimentally confirmed later when the positron 


was discovered in cosmic radiation in 1932. 


1.3. The role of symmetries in physics 


Noether’s theorem states that every symmetry of the action leads to a conserved quantity. 
For example, translation symmetry leads to energy-momentum conservation. Derivations of 
this theorem may be found in many textbooks, and here we only briefly review the following 
example considering invariance under time-translation: Using the Euler-Lagrange e.o.m. we 
get 


lit ggi 


dla OL, DE (4 SL aL . 


d (OL 
=—(—q). 1.33 
dt (Sa) a 
Hence 
d (ƏL 
— (—ġ- L) =0, 1.34 
dt (Za ) Gien 
and the energy of the system 
OL 
BE=(—q—L ie 
(SEA L (1.35) 


remains constant during the motion. 


10 


Energy-momentum tensor of a scalar field. 

Consider a space-time displacement x” — x“ + ĝa, where ĝa is independent from z. 
By computing the change in the Lagrangian £ of a scalar field ¢ theory corresponding 
to this displacement, one easily computes the equations expressing conservation of linear 
momentum and energy: Since £ does not explicitly depend on x, we have 


OL OL 


ôL = m + 5,0) 


Using the fact that [d,0,] = 0 and the Euler-Lagrange equations of motion (1.8), we can 
rewrite this as 


6(Ou@) , do = (0,@)6a” . (1.36) 


OL 
ôL = 0, | == Ov | ba”. 1.37 
» (wan) = 
Alternatively, we may derive 6£ directly as 
OL 
cece 
ôL = azn t : (1.38) 


Since the ĝa” are arbitrary, it hence follows that 


OL 
ð, Tt =0, TY = | ——0,¢ — EL |, 1.39 
MO vy U E (o) v ) ( ) 
where TU is called the energy-momentum tensor. The component Te = (A466 — £) corre- 
sponds to the energy density of the field — cf. Eqn. (1.35). Observe furthermore, that the 


0 component of Eqn. (1.39), 


0 G =F 
at) + VTp = 0, (1.40) 


expresses local conservation of energy, and using the divergence theorem integration over all 
space leads to 


o 
a [eons =0, (1.41) 
assuming the field vanishes at large distances. Similarly, one derives 


qa Fa peor. (1.42) 
expressing conservation of momentum P;. 

Finally, we should mention, that symmetries which are broken at the quantum level lead 
to so-called anomalies. (For example, cancellations of anomalies leads to requirement that 
number of quark and lepton generations match.) 


11 


2. Path integral methods 


2.1. Path integral formulation in quantum mechanics 


We start by recapitulating some basic formulae from quantum mechanics, for simplicity 
consider a one dimensional system, i.e. with one pair of conjugate operators X and P 
satisfying |X, P] = i. The eigenstates of the coordinate operator satisfy X|x) = z|x)} and 
define an orthonormal basis, i.e. 


(x|2’) = (x — 2"), Tiada =1. (2.1) 
Similarly, the eigenstates of the momentum operator satisfy P|p) = p|p) and 


(p|p') =d(p—p'), Jo lp)(p| =1. (2.2) 


The transformation operator between the two basis is given by 


cE = (x|p)*, (2.3) 


1 
(p|2) = T 


and hence the Fourier transform of functions is defined as 
1 l 
Fa) = elh = weleh = em | wto) (2.4) 


The coordinate states in the Heisenberg picture are related to the Schrödinger states as 
|z, t) 7 = e€ |x). Orthonormality and completeness hold for the Heisenberg states only for 
equal times, as can be easily checked. 

In general, when promoting a classical system to a quantum system, implying that ob- 
servables are promoted to operators, it is not clear what has to be done with products of 
non-commuting operators: an ordering prescription, such as ‘normal ordering’ and ‘Weyl 
ordering’ is required. Here we shall consider the latter, also known as ‘symmetric ordering’. 

We may now calculate the transition amplitude U (tp, Te: ti, £i) =H (vp, ty | vi, ti) for 
Lr > ti in the Heisenberg picture. Dividing the time interval between the initial and the 
final time into N equal segments of infinitesimal length < and introducing complete sets of 
coordinate basis states for every intermediate time point, we obtain 


U(ts, Tr: tes £i) = lim dzı -dE N—1 HPE, LIEN 1, LN 1) H 
e—0 


N->0oo 


X H(EN—1,tN—1|EN—2,tN—2) re HLP, HPS, ti) H (2.5) 


with 


5 s TL LPK 
DL elk (tp—tp-1)—ieH (a Pk) 


HlZk,tk|k-1;tk-1)H = (rle "TUH = f E ; 
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using Weyl ordering. Hence, identifying x9 = x; and xy = xf, we arrive at 


dp, ---dpn ee Erai (pr LK THK pk) ) 


U is, Tr: Laa) = lim liee aN 


N= 


(2.7) 


When specializing to the class of Hamiltonians which are quadratic in the momentum vari- 
2 

ables, e.g. H(x,p) = Æ + V(x), the momentum integrals above become Gaussian, i.e. 

particularly simple: 


2 
= | . zarg)? 
[E ($ E = = 2 ( E ) ; (2.8) 
T TE 


hence leading to Feynman’s path integral for the transition amplitude in quantum mechanics 


i zkr)? ktk 
ULL, tf; ti xi) = lim =)” fo . day eo (4 (===) BL s) 


e>0 \2rie 
N-0o 
s tf M 9 
=N | Deexp E dt (Fe -Vv(a)) ) 
= N [De eer. (2.9) 


where N is a constant and S{z] is the action. Here, the end points are held fixed and only 


the intermediate points are integrated over the entire space, which means that the transition 


amplitude is the sum over all paths connecting the two points weighted by the factor e!$!*), 


The dominant contribution is the classical one, namely arising from paths which extremize 
the phase factor, i.e. the ones satisfying aoe 


Exercise 2 Show that for a free particle, the transition amplitude (2.9)) computes to 


m : 
U(t tx) = — Sle] 
(tp wp tea a a ees i 
L U(t, Tf; ti, £i) =ô LT = Ti) S (2.10) 


where T denotes the classical trajectory, and that it solves the Schrodinger equation! 


Exercise 3 Repeat the calculation for the harmonic oscillator! 
2.2. Generating Functional 
Lets us evaluate the matrix element of a product of operators of the form: 


HLP E, tt Xa (t1) XH (ta) |v, ti) g = | dodana arle, H HLT, tı | £2, t2) H 
X H (£2, Hal £i, ti) H 


= Ua (2.11) 
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assuming tı > t2 on the l.h.s. Since z(tı) and x(t2) commute, the path integral on the r.h.s. 
automatically incorporates time ordering. A similar derivation can be done for time ordered 
correlation functions between physical states (cf. B], p.61). 

By introducing the modified action 


tf 
Sle, J] = S[a] + [uss , 


(zf, Lr | Ti, tJ = N f Dz éslre] = (TF; teTeX OIO a, ti) ; (2.12) 
we my rewrite the result above as 


l 1 52 (xp, ty | xi, ti) 
isle] ee 
N | Dev(t)a(ts)¢ 2 57 (hd) = (2.13) 


which justifies the name “generating functional” for (a ¢, ty | £i ti) 7. 


2.3. Path integral for Grassmann variables 


Let 0; with i = 1,2,...,n be a set of Grassmann variables satisfying {0;,0;} = 0. This 
implies that if f(@) is a function of just one Grassmann variable it has the Taylor expansion 


f(0) =a +00 (2.14) 


since 0? = 0. 
A left derivative for Grassmann variables yields 


o 
zg, (19%) = bij 9k = ikl; s (2.15) 


Notice also that E = 0, i.e. are nilpotent as well. Hence, E i 9; } = bij. 

The notion of integration can also be generalized to Grassmann variables. Considering 
that the integral of a total derivative should vanish if we ignore surface terms and that 
a definite integral, being independent of the variable, must give zero upon differentiation. 
Since differentiation with respect to a Grassmann variable is nilpotent, it satisfies these 
properties and integration with respect to Grassmann variables can be naturally identified 
with differentiation, i.e. 


[oro m on) | (2.16) 
Therefore 
[=o [0 =1, (2.17) 
md 


f dof (0) =a f 49 PI Ja) (2.18) 
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for 0’ = aĝ. Generalizing to several Grassmann variables we find 
fa ...dOn f (0i) = (det ai;) fa dolf (a 8) 


for 0; = aij; and (det a;;) A 0. Furthermore, a delta function can be defined as 


6(0) =0 = —i f dis , 


where ¢ is another Grassmann variable, since 


[950 =f, e610) 0) = (0). 
Finally, a Gaussian integral for Grassmann variables is given by 


*/ pf + ML 
J] «978; (FF Mis LGT i+ 07 C1) — det Mi; Los 0. ei bitce] M,j cj 
1.1 


1. 
= const. det i: e% Mij © i 


assuming ĝ0ž and ĝ; as well as c¥ and c; are independent Grassmann variables. 


Exercise 4 Consider the following action for the fermionic harmonic oscillator: 


Sto. Panl = Sd, 8] + f iat + d 
om H 2 GF = 
Stu.) = fat (2 - by) -= 0 ; 
and compute the free fermion path integral! 


2.4. Path integral for a relativistic scalar field theory 


Let us consider the action of ¢4 theory 
d = fatec(o.9,0), 
= 49,60") — ge Apt, 
and introduce appropriate source through the couplings 
S[d, J] = Sl + I d 


leading to the e.o.m. 
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(2.19) 


(2.20) 


(2.21) 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


Furthermore, considering the limits t; -— —oo and ts + +00, the vacuum to vacuum 
transition amplitude is given by 


Z[J] = (0|0)7 = N f Does (2.27) 
For A = 0 we easily derive 


Zul =N f Doele] 


=N (det( +m?) as [dtad*a! J(x)G(a—a')J(a') 


— Zo[0]e72 TdT J(x2)G(x—a') I(x’) (2.28) 


where the field ¢ was integrated out by shift and quadratic completion, Solo] = S|] | = 
and G(x — x’) is the Feynman propagator resp. the Green function we encountered earlier 
in Eqn. (1.20). It is therefore straightforward to see that the normalized time ordered 
correlation function 


© 1 PZT] 
a20 PZJ] od (x5. (2) 


=iG(x— x’). (2.29) 
J=0 
The path integral Z[J] of the full theory with 4 0 cannot be evaluated in a closed form, 
but for weak coupling A one can compute it perturbatively. By making the replacement 
LT) > a in the interaction part of the action, we can derive the following power series 
in the coupling A: 


ZA = (<* J (ah) ) N I Ddu = ts J (eh) ) Zol J] 


(OlT 6(x)d(2’)|0) 


iA fa & lf ANT Pa 0 a, 6 
=j y PTA EPT dr = + ...| Z 
t= J TEJE) 2 ( 0 J h | tapay a 9 
iv 07 i PAA g4 
=Z Li 4 al es Jd xid r2J(x1)G(x1—x2)J (x2) l 2, 
K h T f "J" Je ? (2.30) 


Let us compute the 2-point correlation function to linear order in À. For this we need 


Z| J] = Zo[0]e~2 fd4x1d*a2J(x1)G(x1—2x2) J (x2) 


1+ P G r(0)Ge(0) fe 


+ Zaro) la, — z3)J(x3)Gp (z — 24) J (24) 


7 £ d'xd‘x3_¢G(x — 13) J(x3)G(x — z4) J(£4)G(£ — z5)J(x5)G(x — xg) J (z6) 
Ñ 003 (2.31) 


which follows from Eqn. (2.30) by explicit computation while consistently keeping only terms 
linear in A. Notice that 


Z[0] = Zo[0] (1 + 6 r(0)Gr(0) fe) L OLA. (2.32) 
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is a divergent quantity which can be absorbed in the normalization M. Continuing our 
computation of the two-point function to order À, we arrive at 


TT 
"aeaa = GZ TTA ed 1... 
= iGr(x1 — £2) — *Gr(0) [accra — 21)Gr(x — 22) + O(A?). 


(2.33) 


Similarly, we may compute the lowest order quantum correction to the 4-point function, 
which however involves numerous terms. Therefore, a systematic procedure keeping track 
of the perturbative expansion is helpful. Such a procedure is given by the so-called Feynman 
rules: 


Ly T2 iGp (z1 — Ta) 
Tı J T4 l TSI 

E V (z1, £2, £3, £4) = EREET (z3) (x4) 
T2 T3 


= —ià\ frase — x1)ðt (x — z2)ð (x — z3)ð (x — 24). 
(2.34) 


With these rules, supplemented by the rules that we must integrate over intermediate points 
in a graph and that if the internal part of a diagram has a symmetry one must divide by 
that symmetry factor] we find for the simplest 2-point diagram 


1 . i c 
i A S [ences — yi )iG(y2 = x2)iG(y3 = ya) V (41, Y2, Y3, y4) 


= -2G(0) [evo ~ WG — 22). (2.35) 


Hence, we see that the result of Eqn. (2.33) can be written diagrammatically as a simple 
sum: 


(UIT 9(x1)G(x2)|0) = + | | +... 
Tı T2 


Tn order to compute the correct factor for a certain type of Feynman graph, one first counts the number of 
possibilities to connect the external lines to the vertices of a “pre-graph”. That number is then multiplied 
by the number of possibilities to connect the vertices among each other, and the result is divided by 4! for 
every vertex and by n! if n identical vertices are present. The latter n!-factor comes from the expansion 
of the exponential in the path integral. In the present 2-point graph, one arrives at 4.3/4! = 1/2. 
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When computing the first order corrections to the 4-point function, we easily see that 
connected and disconnected graphs appear in the result of 


4 
(0|T(@1)o(w2)o(w3)o(ea)|0) = == Z|] 


D 63 (a1)6I TEA (aa) (2.36) 


J=0 


T1 
: x La x 2 
+ 6 permutations of ia t L >< TOI). 
i.e. the generating functional Z [J] generates both types of graphs. In contrast, the Z,[J] := 
—iln Z[J] generates only connected graphs and is hence called the generating functional of 
the connected graphs. For example, observe that 


1 ZJ] (i ( A- PZ) A 
i SJ (æ1)ðJ (22) ZH) ZI OI (ar) 6T(w2)) | 
= (0|[TØ(z1)$(z2)|0) — (0[(21)|0) (0| (x2)|0) 
(OlT E(z1)$(z2)l0)c . (2.37) 


In the present case, however, (0|¢(x1)|0) = 0. But for e.g. the 4-point functions there is a 
difference, and one can convince oneself that 


1d) l a 
j4 6] (x1)ôJ (x2)5I (£3)ôJ (£4) = (UIT hL 110 2) P( 3)@( 4)l0)c = A S (2.38) 


Exercise 5 Explicitly compute the 4-point function up to order A and show that Ze generates 
only connected graphs while Z generates both connected and disconnected ones! 


At order \? one discovers, that the connected diagrams generated by Ze still contain two 
types of diagrams: the one particle reducible ones, which can be reduced to two connected 
diagrams by cutting an internal line, and the more fundamental one particle irreducible 
(1PI) graphs. 


2.5. Effective action 


Let us define the “classical field” $,(x) a$} 


belz) = W160) = asst = A, (2.39) 


Note, that (x) is indeed a functional of the source J and in the present case of ¢* theory vanishes for 
J = 0 — in general it becomes a constant. 
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The name “classical” is motivated by the following observation: Assuming that the measure 
of the path integral does not change under a redefinition of the field variable, we have 


SZ J) =N [pei l d'aso) los [J] <9, (2.40) 


since Z is independent of ¢, and it follows that the Euler-Lagrange equations hold as an 
expectation value equation (Ehrenfest’s theorem): 


ôS |, J] 


OT IW: =0. (2.41) 
Considering the generic form ER) = F(¢(2)) — J(x) we find 
-v foo Bea aston =0, or 
L (a E J) z=, 6 
F (F — sno) — J(z) =0. (2.42) 


Restoring A (which have set to 1) in this equation and using Eqn. (/2.39), we get 


F u = ra — T (2.43) 


which in the classical limit A — 0 reduces to the form of the classical Euler-Lagrange 
equations. 

The relation indicates that the variables J(x) and LT] are in some sense con- 
jugate variables. This motivates the definition of a new functional through the Legendre 
transformation 


Tio] = Z,|4| = f d'el). (2.44) 


Clearly, oes = —J(x) which reminds us of se = —J(x), and therefore I'[¢,] is known 
as the “effective action functional”. Note that 


ja _ 8 ZelJ] FT [dc] 


miaa ln T (2.45) 


which follows from = = —ô(x — y) using the chain rule and the definition of Qe. 


Introducing the abbreviations 


= Spela) ETT fE T SIE) 6 En) (2.46) 


one can show that ac: Ñ ‚i.e. P™ evaluated at [J = 0], which for ¢*-theory is 0, is the 
1PI n-point vertex function. Hence, I'[¢,] is also known as the 1PI generating functional. 
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For example, since Vial = —G?) Eqn. (2.45) tells us that PG, is the inverse of the 
full propagator at every order of perturbation theory. Hence 


1 
GO) = — ~ = Gp+GrUGrp+GprlGruGrt..., (2.47) 
Gr -È 
using Eqn. (2.45) in a short-hand notation in the second step. Furthermore, it follows from 
varying Eqn. (2.45) with respect to J (w) that 


fe 8? ZelJ] STe] _ L Zd) FT [dc] 8? Zel] 
SJ (w)dI(y)bI (z) ôpe(z)ðpe(1) J (y)ðJ (2) dbc(z)dbc(x)dbc() SE 
2.48 
and hence using once more one finds 
"ZJ = Ln 07Z-1JT 07T ô? Ze[J] ST [be] 
J (w)ðJ (y)ðJ (z) ôJ (x) (x) 8J (y)bI (y) J (2)8J (2") ee 
2.49 
or in short hand notation 
Z® = ZØ ZØ ZOTO. (2.50) 


Similar relations can be derived for the 4 or n-point Green functions by varying Eqn. (2.45), 
i.e. for the former one has in graphical notation: 


C K (peat 


2.6. S matrix 


Let us define an incoming free field ¢j,(z) := lim assuming an adiabatic vanishing of the 
xLQ—— Co 


source j at zo — +00. Similarly, an outgoing free field is then @out(x) := lim . The idea 
TOK -too 


is, that a free incoming field propagates in time, a source is adiabatically switched on, the 
field interacts, and the source is adiabatically turned off leaving a once more free outgoing 
field. This construction, of course takes place in a given Hilbert space, and one may try to 
find the unitary operator S which connects the in- and out-fields (and states): 


dout(2) = B". lout) = Sin) . (2.51) 
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In general, one is of course interested in the situation where an initial configuration of 
particles œ ends up as a final configuration 8 after a scattering process. The scattering 
amplitude for this is denoted 


Sap = outl 8 | @)in = in(8|S|@) in : (2.52) 


and its square measures the probability for the process in question. If for example |a)in 
consists of two scalar particles with momenta pı and pg, then |@) in = al, (pi)at, (p2)|0) where 
a' are creation operators. Obviously, invariance of the vacuum demands that Soo = 1. 

In Eqn. we have already given a general solution to the classical e.o.m. for a scalar 
field. Using the expressions for advanced and retarded Green functions of we may 


write 
olz) = Pin(z) ale lage. HIE = Pout (T) + l IIE : (2.53) 


and consider the weak asymptotic condition Jim (a|@(x)|b) = (a|@our(x)|b). (Note, that we 
T00 in 


are considering quantized fields ¢ in this section.) 
Now, the inhomogeneity j above is due to the interaction part, i.e. should be replaced by 


jl) > = Ky¢(z'), (2.54) 


where the last step follows from the e.o.m and K = (O+ m7?) is the Klein-Gordon operator. 
Next, we define the functional 


I[J] = Teil dad (x)o(x) (2.55) 


whose vacuum expectation value is the generating functional Z[J], i.e. (O|Z[J]|0) = Z[J]. It 


then follows from (2.53), and S—1 = St that 


Pout T — din = i fate (Gret(x, 2’) — Gaay(x, 2")) Key Si ; and 
lm); SI[J]] = i fate (Gret(x, 2’) — Gaay(x, 2’) Key (2.56) 


Now one can show that upon promoting #* in Eqn. to creation/annihilation oper- 
ators, that the expression Gadv — Gret coincides with the commutator of free scalar fields 
dn LT), Gin(x’)]. Employing furthermore the Baker-Campbell-Hausdorff formula in the form 
[A, e?] = [A, B] e? for the case where [A, B] is a c-number, the solution 


Sr = exp ( olarsa lF 


where F[J] is some functional of J, suggests itself. Taking the vacuum expectation value 
thereof and considering (0|.S = (0| in the absence of external fields, we find (for normal 
ordered ST): 


(O|SI[J]|0) = F[J], and (0|SI[J]|0) = (O\Z[J]|0) = Z[J]. (2.57) 
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Hence, we arrive at the so-called reduction formula 


g= exp( [ots én) Koga | ZU) 


where “:” denotes normal ordering. The above formula means, that for a given Feynman 
graph, one has to amputate all external legs and replace them by the incoming free field wave 
functions in. (The operator K converts external propagators into 6-functions.) Hence, the 
n-particle S-matrix element is 


(2.58) 
J=0 


H = | 1 6G) eC (2.59) 
i 
where G(21,...,2n) denotes the n-particle Green function. 


2.7. Euclidean path integrals 


By analytic continuation of all integrals to imaginary time in the complex t-plane, i.e. 
t => t’ := —i7 where 7 € R, one arrives at Euclidean path integrals. This “Wick-rotation” 
can be useful to simplify certain computations, taking into account the famous Oster- 
walder—Schrader theorem which basically states that if a quantum field theory exists in 
Euclidean space, it also exists in the Wick-rotated Minkowski space-time. Of course, one 
must take care not to cross any poles when doing a Wick-rotation. For example, consider 


3(k°) 


Figure 2.1.: Wick-rotation 


the following propagator (of the harmonic oscillator): 


1 1 1 
Gr(k) = ln ~=, 2.60 
r(f) paar A RTU iek —wt+ie ( ) 
which has poles at k = —w+ie and k = w—ie. Hence, for the analytic continuation from Rk° 


to Sk° one must rotate anticlockwise, i.e. k? > k” = ix, as indicated in Figure|2.1} Since 
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we can represent k? —> ig, it follows that in the complex t-plane the consistent rotation 
will be t — t = —ir as stated above. 

Euclidean path integrals can in fact be interpreted in a statistical context. Let us make 
this clearer by reviewing some facts of statistical ensembles: For a given ensemble (e.g. 
a system of harmonic oscillators), the value of any observable quantity averaged over the 
entire ensemble will take the form 


(A) = X pn(n|Aln) = X pada ; (2.61) 


where pn is the probability of finding a system in the ensemble to be in an energy eigenstate 
|n). In fact, there are two kinds of averaging involved here: the average in a quantum state 
(the expectation value) and the averaging with respect to the probability distribution of 
systems in the ensemble. As pn is a probability, it has to satisfy the conditions 


0< pm <1, X i=l (2.62) 


n 


If we consider a thermodynamic ensemble interacting with a large heat bath, and assume we 
have waited long enough to achieve thermal equilibrium, then the probability distribution 
is given by the Maxwell-Boltzmann distribution 


1 1 
—BEn = 2. 


where E is the energy of the nth quantum state, kg is Boltzmann’s constant, and T 
the temperature of the system. Using the conditions for ppn, one easily determines the 
normalization Z: 


LS anl ae L and (2.64) 


n 


In fact, Z (8) is known as the partition function of the system and plays the most fundamental 
role in deriving the thermodynamic properties of the system. 
For such a thermodynamical ensemble, the thermodynamic average according to (2.61) is 
given by 
1 -6H Tr (e941 A) 
Ls = sm > (nle FP" Aln) = aoe (2.65) 


In particular, the average energy computes to 


usetan e (2.66) 


Z(6) op 07 
and entropy is given by 


Unna = S = Ñ "pa ln Ba = X pn (BEn +n Z(8)) = BU + n Z(B) 


— -i (5216) (2.67) 
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The free energy can hence be written as 
FLH =U —kgTS 
1 


and hence the partition function takes a particularly simple form when expressed in terms 
of free energy: 


Z(B) =e PF) | (2.69) 


To sum up, a thermodynamical statistical ensemble can be described by a Euclidean path 
integral where the Euclidean time interval plays the role of temperature and the partition 
function is identified with the transition amplitude. In general, one uses Euclidean path 
integrals to study phase transitions of quantum systems — an example from high energy 
physics would be a quark-gluon plasma. 


2.8. Path integrals for gauge theories 


Our starting point is the generating functional 
ZL = eel = N I DAES lAn], (2.70) 
where 
Sz[Ayl = f d'z (irer + jua) 


= : [ae [irao — x’) A(z") + Tan. 
OM (a 01 = (Cnt? — 846") 5a — x’), (2.71) 


neglecting surface terms. The path integral over the gauge fields A, would lead to det O. 
Unfortunately, this integral does not exist since O is a projection operator (cf. Section|1.2.2). 
In fact, the source of this difficulty is the invariance of the above action under a gauge 
transformation 


A, > A) = U ALUT +iU* (ð U), 
U (a) = e0) , (2.72) 
All Al) that can be obtained from a certain A, by making a gauge transformation with 
arbitrary a(x) are said to lie on an “orbit” in group space (in this case the Abelian U(1) 


group). Since the action is invariant on such orbits, the generating functional Z[J„] is 
proportional to the “volume” of the orbits denoted by f TT da() (resp. the group invariant 


T 
Haar measure || dU (x) in the non-Abelian case). This infinite factor needs to be extracted 


T 
before doing any calculations, and this is best done by the method of Faddeev and Popov: 
The idea is to “fix” the gauge freedom and to integrate over each orbit only once, i.e. one 
chooses a hypersurface defined by a “gauge condition” F(A,) = 0, which intersects each 
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gauge orbit only once. Even if A, does not satisfy this condition, one can find a gauge 


transformed Ale a) which does, i.e. F F(A) = = 0 then has a unique solutio E for a(x). 
Now the trick due to Faddeev and Popov to extract out the infinite gauge volume factor 


relies on the insertion of 
App|A n J Tlc) ö (F(A Alo PRT = =I; (2.73) 


into the path integral. Note that Arp[A,] is a gauge invariant quantity because its inverse 
is. (The measure in the group space is invariant under a gauge transformation.) We make an 


inverse gauge transformation A,, —> AL. absorb the gauge volume into the normalization 
N and arrive at 


=N f DA, Arp[A,]6 (F(A,,)) RA (2.74) 


In order to determine Arp|A,], notice first that 


= a= J Tee) (Feat a= JIar (r P(A) («er 5) 


ôa 
= det : 2.75 
Thus, we find the Faddeev-Popov P 
5F (AM 
esta aa Sed (2.76) 
a a(x)=0 


which can be thought of as the Jacobian that goes with a particular gauge choice. 

In a further step, we generalize the above to F(A,(a)) = f(x), where f(x) is independent 
of A,. Since physical quantities are independent of f, we may multiply the generating 
functional by a weight factor and integrate over all f (x) to arrive at 


=N [RABLE Hl RA Se 
-N 1 DA, Apples A A Ste Pane) (2.77) 


where € is known as the gauge fixing parameter. In fact, this new term in the action can 
equivalently be seen as the result of integrating out a new multiplier field b. Namely, observe 
that 


N J Dbe laL) — p'e R SAE A. (2.78) 
Finally, we can also write the Faddeev-Popov determinant in a path integral form by intro- 
ducing unphysical “ghost” fields c and €: 


P(A) (x) 


—i fdtx dar) (ae) GLC) 
Arp|A = | Pence a=0  , (2.79) 


3One should mention at this point, that such a hypersurface cannot be found in the non-Abelian case in 
the sense that any hypersurface will intersect gauge orbits more than once. This is commonly referred to 
as the Gribov ambiguity. 
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Thus the complete action (without sources) for this choice of gauge condition reads 
1 
S[A,, 6, é, c] = [ae (ABP + (O,A")b + Sy + d.ra : (2.80) 


dropping surface terms once more. In the present model, the ghosts do not interact with 
the gauge fields and can hence be neglected in explicit loop computations, however, in 
non-Abelian gauge theories this is no longer true. 


Exercise 6 Show that the gauge field propagator in momentum space derived from the ac- 


tion (2.80) is given by 
=i PuP 
Gyv(p) =o (w = (1 Ñ ek - (2.81) 
p D 
Why is it impossible for the auxiliary field b to appear in Feynman graphs? 


BRST invariance. 

Since we have actually “merely” inserted a factor of unity, Eqn. (2.73), into the path 
integral, the physical content of the theory has not changed. On the other hand, the action 
Eqn. no longer exhibits the gauge symmetry due to the gauge fixing terms. However, 
this action exhibits a new global nilpotent fermionic symmetry, which in some sense “re- 
members” the gauge symmetry of the original theory. It is called the BRST symmetryļ4, and 
in our case the according symmetry transformations which leave the action above invariant 
are given by 


sA, = Ouc, sc=0, 
SC <D. sb= 0, 
a’ =0. (2.82) 


If one repeats the same Faddeev-Popov procedure for non-Abelian SU(N) gauge fields 
where 


Fup = Oy Ap — Op Ay — ig (Ap, Av] = FLT", 
Fi, =0,A% — 0, AG + gf PAPAS, 
| =i 
Sinv = I dir= TF F = f d's- FaFa, (2.83) 


f% are the antisymmetric structure constants, and T" are the generators of the SU(N) 
group, one arrives at 


3 


1 a Qa, uV a a a =a a 
Shed = fale (PRL Pm + MAB ESOP + De) , 


De = 6, + gf AS. (2.84) 


‘This symmetry was discovered in 1975 by C. Becchi, A. Rouet and R. Stora, and independently by I. V. 
Tyutin in the same year. 
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This action is invariant under the BRST transformations 


sA; = (D,zc)*, so = See 
sc’ =b", sbt =0, 
2? =0. (2.85) 


Note, that as before the BRST transformation of the gauge field has the same form as 
the according infinitesimal gauge transformation, but with the ghost c(x) replaced by some 
(bosonic) function a(x), i.e. 


Anla) > UHA HI + UH) = Aula) + Dala), 


U(x) = eTA OT" (2.86) 


Furthermore, in the non-Abelian case the field tensor transforms covariantly under an in- 
finitesimal gauge transformation, i.e. as 


Fy > UFU & Fuy — ig [Fuv, a(2)] , (2.87) 


(in contrast to the Abelian case where it is invariant). However, the trace of two field 
tensors TYE... Fp (and hence the invariant part of the action) is gauge invariant due to 
cyclic invariance of the trace. 

As a final comment to this section, we mention that the physical space of the gauge theory 
can be identified as satisfying 


Qprst|phys) = 0, (2.88) 


where Qprgrt is the charge constructed from the Noether current of the BRST symmetry. 
Additionally, physical states are free of ghosts, as should already be clear by how the ghosts 
were introduced into the theory. Since Uns = 0, our Hilbert space actually decomposes 
into three subspaces: one where Qprst|w) 4 0, one where states may be written as |Y) = 
Qprst|v’) and hence are annihilated by Qprsr, and the physical one above but where 
states cannot be written as in the second subspace. Furthermore, the S-matrix of the 
physical subspace must be unitary so as not to mix physical and unphysical asymptotic 
states. This is of course the case, as can be shown. 

In the Abelian case, condition implies the Gupta-Bleuler condition 0“ A,,|phys) = 0 
for BC fixing we considered above. One may also verify through explicit computation 
that s = 0 using Eqn. and the fact that the vacuum belongs to the physical 
Hilbert space. 


Ward identities. 

BRST invariance of the theory leads to various relations between Feynman graphs (and 
hence scattering amplitudes) called Ward identities. These play an important role when 
it comes to proving renormalizability, and additionally provide various consistency checks 
when doing explicit computations. 

In the derivation above, we have introduced several new fields into the gauge theory for 
which we need to introduce sources as well. we consider the non-Abelian case (of which the 
Abelian theory is a simpler special case), and write 


SAR €, c] = S[A,b,é, c] + ae (JOM AS + 5204 e — En?) . (2.89) 
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A BRST transformation of this action inside the path integral leads to 
= = 4 “A, [L a I pabc-a bc _ paa 
sS7[A,b,¢,c] = | dtz (j (Dpue) rS NEE — bne), 
sZ =0=N f DADbDeDei (sSJ|A,b, Z, c]) 7/4641 
= dł Fah DH a -a g abc b 6% _ „a hn 2.90 
x LPD GI + ES FO ere) — PEH) (2.90) 
If we introduce two further external sources for the non-linear BRST transformations 


Sext = [ae (Ko D0)" = E eee) l 


SJ, KlA,b, c, c] := SJ[A,b, €, c] + Sext» (2.91) 
we can rewrite (2.90) as 
. ZJ, K]  _ ôZe|J, K] 6Z,|J, K] 
4 a, mimo: 7 es L H = U. 2.92 
je T (: (x) 5K H(z) n (x) AKT) n (x) UE U ( 9 ) 


From this “master equation” one can derive all identities relating the connected Green 
functions of the theory. A Legendre transformation allows to rewrite the master equation 
in terms of the effective action, i.e. with 


VAL. ¢,c, K] = ZJ, K] — [as (GOP AG + jD + HC? — en’) , (2.93) 


(where we dropped the subscripts c of the classical fields — cf. Section |2.5)) it follows that 


= 4 or or AT or gp. E 
k (ne dK G(x) N dcn) K(x) + Tea)” ) 0. (2.94) 


In this form, the master equation is commonly referred to as the Slavnov-Taylor identity, and 
it allows to derive all relation between the 1PI vertices resulting from BRST invariance of 
the theory. Note, that at tree level (i.e. in the classical approximation) where T| ġe] = S|}, 
the Slavnoy-Taylor identity is just another way of writing sS[¢] = 0. 
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3. Renormalization 


3.1. Regularization and power counting 


In general, perturbative computations in quantum field theories involve distributions multi- 
plied with one another. These expressions are mathematically ill-defined and hence usually 
lead to divergences. In order to handle such expressions, the common strategy is to regu- 
larize the expressions in order to be able to split off the divergent parts which then may be 
absorbed by redefinition of the parameters of the Lagrangian. The latter process is referred 
to as renormalization. 

Let us illustrate this procedure by considering a scalar quantum field theory, namely ¢4 
theory: From Eqn. we see that the propagator diverges for 7 = 2’, i.e. 

dtp 1 
G(0) = — f (2m)! p? — m? ` (3.1) 

Since there are four powers of p in the numerator and two in the denominator, this integral is 
expected to diverge quadratically at large p (hence the term “ultraviolet/UV divergence” ). 
Similarly, the square of a propagator is a divergent quantity as well: 


N2 — dp ip(e—a! dk o 
ee-e- fam [eRe C 


which diverges logarithmically for large k. 


Power counting. 

The superficial degree of divergence of an arbitrary Feynman graph can be inferred by 
considering the properties of the Feynman rules for propagators and vertices of the model 
in question. For ¢* theory, for example, one has one propagator behaving like 1/k? for large 
momenta, hence reducing the degree of divergence by 2. Furthermore, every 4-dimensional 
loop integral raises the degree of divergence by 4. Hence, 


diah =4L — 21, (3.3) 


where L denotes the number of loops and I denotes the number of internal lines (i.e. prop- 
agators). Since there are J internal momenta as well as momentum conservation at each 
vertex and finally also overall momentum conservation, the number of independent mo- 
menta, which is L, is given by the relation 


pa =i (3.4) 


where V denotes the number of vertices in the graph. Finally, one needs relations between 
the number of vertices and the number of legs. External legs denoted by E count once, 
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whereas internal legs J count twice as they are always connected to two vertices (or two legs 
at one vertex). Thus, one finds the following relation: 


4V=E+21. (3.5) 
Putting all pieces together, we hence find 


d(y) =4+ 21 — 4V 
=4-E. (3.6) 


Notice, that the simultaneous elimination of J and V was only possible in 4-dimensional 
space-time. Furthermore, the mass dimension of the scalar field @ is 1, as can be seen from 
the action Eqn. (2.24). In general, the field dimension will be related to the numerical 
prefactor of E in the according power counting formula. 


Exercise 7 Check this claim by deriving the according power counting for QED! 


The power counting formula for scalar ¢4 theory Eqn. tells us, that any correction to 
the two-point function can at most diverge quadratically (E = 2 — d(y) = 2), and any cor- 
rection to the four point function can at most diverge logarithmically (E = 4 — d(y) =0) 
independent of the number of loops. The fact that only a finite number of n-point functions 
(namely two) exhibit divergences is an important requirement for renormalizability of the 
model, since we only have a limited number of parameters we may redefine to absorb these 
divergences. In order to do so, non-existent divergent integrals must first be regularized. 

Various regularizations schemes are known, and we start by reviewing Pauli-Villars reg- 
ularization applied to the example above and then compare to dimensional regularization. 
Both schemes preserve translation and rotation invariance of the theory. 


Pauli-Villars regularization. 
This regularization scheme was introduced by in 1949 by W. Pauli and F. Villars and is 
based on substituting the Feynman propagator G p(x — YT. m) by an expression of the form 


GOLE 01 = X iGr(x —2',M;), 
i=0 
co = li; Mo=™, (3.7) 


supplemented by the conditions 


yami, (3.8a) 


y aM =0. (3.8b) 


The first condition (3.8a}) cancels the strongest singularities of any Feynman graph, and 
the second one (3.8b) then cancels remaining singularities. This means, that if a given 


theory is only logarithmically divergent, condition (3.8al) suffices, and only one auxiliary 
mass Mı = M with cı = —1 is necessary. 
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Let us illustrate how the first condition removes the quadratic singularity in G(0) of 7 


theory: 
4 
-j d'p clle) ee 
(27)4 p-m? TE 


=- | ant Gm MA’ ee) 


GEE- 2") 


which for x = 2’ is only logarithmically singular instead of quadratically. Obviously, all 
divergences are recovered in the limits M; —> co, i> 0. 

Let us compute an exemplary graph in $4 theory: 

Allin all, Pauli-Villars regularization works fine in scalar field theories and also in QED, 
but it breaks gauge invariance when considering non-Abelian gauge theories such as QCD. 
Therefore, dimensional regularization is usually the best choice in those models, as it always 
preserves gauge invariance. We shall review its properties in the following. 


Dimensional regularization. 

This regularization scheme was introduced in 1971 by G. ’t Hooft and M. J. G. Veltman 
and is based on the observation that divergent Feynman integrals would become convergent 
when computed in a smaller space dimension. Therefore, computations are done in arbitrary 
d-dimensional space and the result of the integration can then be analytically continued 
to real or even complex values of d. The original ultraviolet divergences then manifest 
themselves as poles at d = 4, typically in the form of Gamma functions. 

We illustrate this procedure by repeating our previous example in this scheme: Since 
we are now considering 4 — < dimensional space-time, the canonical dimension of the cou- 
pling A would change as well, unless we introduce a new parameter ÀA —> uà. Hence, 
using Eqn. in Appendix[A.3 the one-loop correction to the two-point function in 


momentum space is given by 


l d‘-€p 1 idm? (Arp? SZ 
= p rre C T(-l 2 
wr fa a (-1+6/2) 


2 2r)4—E p? — m2 —m? 
id 2 
N — + finite. (3.10) 


Observe, that the above integral becomes zero in the limit m — 0. This is a typical feature 
of dimensional regularization. 


3.2. Renormalization and renormalization group 


Mass renormalization. 

The result we have just obtained in (3.10), is essentially what we denoted earlier as —iX 
— seeEqn. (2.47). Defining the “physical” (or renormalized) mass mphys by the pole of the 
complete propagator 


(3.11) 
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gives on comparison 
2 wg = + E=m? (1 a (3.12) 
Mohys m‘ + m , . 
which to order A is equivalent to 
ae (3.13) 
pays 16r?e) ’ 


where the renormalized mass is given by m 
2 


ao = —I(?)(0) and is taken to be finite (i.e. the 
bare mass m* is infinite and compensates the 1-loop correction). 

An alternative point of view which is quite common is to consider counter terms in the 
Lagrangian and treat those as an interaction giving rise to additional Feynman rules. In the 
present case such a counter term would read 6£L; = —16m?¢? and its effect is that we now 
have SIL) = k? — m? since the divergence is cancelled by the additional “interaction”. 
Hence, in this picture m? is a finite quantity (in contrast to before). The introduction of 
this counter term is equivalent to multiplying m by a renormalization factor Zm. 

The coupling 4 is renormalized in a similar way by computing the 1-loop correction to 
the four point function, i.e. Aren. = re). 


Exercise 8 Compute the 1-loop correction to the four-point Green function, showing that 
in dimensional regularisation it results to 


3A 


Dip) = iue [1 -2 
DY (pi) = LAU (1 ic 


) + finite. (3.14) 
Renormalization of the wave function ¢ is not necessary at 1-loop level in this model, but 
it appears at two-loop level. We call the according renormalization factors Z) and Z¢. 


Renormalization group. 
We have the following relation between the n-particle vertex function and its renormalized 
counter part: 


T (pim, A) = Z3 (ape (pi, mr, Ar, u) - (3.15) 
It the follows that 
o 
— pm = 3.16 
Ma, 0, (3.16) 
and hence 
o o OA, O Om, O 
—np In \/Z5 + p= L S F =0. 3.17 
( "Hap ale oF E Dy I Aau i) en) 
Defining the quantities 
Om, 
r= — | Z. ; m Àr = , 
IAr) T ™Mm(Ar) E 
Or, 
BA) = 3.18 
GO = 0 (3.18) 
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we may write this equation as 


o ð o 
(2 + BOr) ay = ny(Ar) F Than LM im) ro”) =0. (3.19) 
It is called the renormalization group equation and expresses the invariance of the renor- 
malized I) under a change of regularization parameter L- 
Let us now derive a similar equation expressing the invariance of r) under a change of 
scale, i.e. p > tp, m > tm, u > tu. Since In] has mass dimension D = 4- n +€ (3 — 1), 


r (tp, A,m, u) = tPr™ (p, à, m/t, w/t), 


a a a 
a rae (n) 
(+5 +m tH p)r 0, (3.20) 
and using (3.19), 
Aa eee j= 4. Dd T™ (tp,,m,p) = 0 (3.21) 
OA G Ym D.T, U) 5U, - 


where we have omitted the subscripts r. This equation tells us that a change in t may be 
compensated by a change in m and X. Hence, we expect the form 


(tp, Am, u) = FOL (p, A(t), mt), 10 (3.22) 
and differentiating this with respect to t reveals that +20 = B(A). Therefore, A(t) is 
called a running coupling constant. The zeros of the -function are called fixed points and 


depending on the quantum theory several of these may be present. 


Exercise 9 Compute the 1-loop 8-function for * theory and show that it results to 


tes, ON 3X 
B(A) = lim PS Sre. >0. (3.23) 
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A. Supplemental Material 


A.1. Dirac formalism 


In this section, we will give a short introduction to the quantization-formalism introduced 
in 1964 by P.A.M. Dirac which enables to handle constrained Hamiltonian systems. To this 
end we will mainly follow reference [6]. 


A.1.1. Hamilton systems with constraints 


We start with an action for classical mechanics: 
S= J Han dn)at, (A.1) 


where L is the Lagrangian and the da denote the time-derivatives of the generalized coordi- 
nates gn. Variation of this action leads to the Euler-Lagrange equations 


d (OL OL 
— | — | = —, A.2 
dt (sn) diu a 
and using the chain rule for the left hand side of this equation leads to 
- PL dL OL 
Gn \ > 5. h «dC P a Un (A.3) 
On! OGn odn OGn! On 


da: can only be determined from this equation if the determinant 


07 L 
A.A 
det C #0, (A.4) 


is unequal zero. In this case, a Legendre transformation leads to the so-called Hamiltonian 
H (qn, Pn) = Pnån (dn, Pn) — Llan, dudu: Pn)), (A.5) 


where the conjugate momenta pn are defined as 


OL 


Dan x (A.6) 


Pn = 


If the pn do not depend on the gp, which is exactly the case when the determinant (A.J 
vanished}, one gets certain relations 


Pm(4 p) =0, (A.7) 


% 2 
1With {A.6) one has =% = 22, 
Ody! Odn ån! 
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out of (A.6), so-called primary constraints. But even if the transformation is singular, one 
can easily show that H depends only on qn and pn: Variation of the right hand side of (A.5) 
yields 


. z OL OL .. ; OL 
ÔH = prdgn + GnOPn — ag = Td = GnOPn — 3 Tin. (A.8) 


where the definition (A.6) was used. Hence, from (A.5) follows 


OH OL OH 
la = dpa — Fae t+ San) bn = 0. A.9 
(4 a) á (= Odn Ñ ( ) 


Note that the variations ôpn and ôqn are not independent from each other because of the 
constraints. Since every function G on the phase space, which vanishes on the subspace 
m = 0, can be written as a linear combination of the constraints (G = gmm), one concludes 
that (A.9) must have the form 


0 0 
tn Soe Spp + tn 5 Sgm =); (A.10) 
(A proof can be found in e.g. reference [6].) Comparing coefficients finally leads to the 
generalized Hamiltonian equations of motion 


OH Obm 
n = D7 + Um ; A.lla 
Opn Opn l ) 
OH Obm 
a =- — Unm : A.11b 
j On On ( ) 
where the equations (A.2)) and were used for the left hand side of (A.11b). 
We will now need the definition of the Poisson bracket 
Of ð Of O 
{f, g}rs = f mae eee (A.12) 


qn pn OPn On’ 


which is antisymmetric and fulfills the Jacobi identity. Let g be an arbitrary function of qn 
and pn. Its time derivative is then given by 


l) Og . Og . 
= e n + Pn. A.13 
Dan. Opn? ( ) 


Using the Hamiltonian equations of motion (A.11) and the definition of the Poisson bracket 


(A.12}) one obtains 
g= {9 A} pp ka Um{9, Pm trs . (A.14) 


Following the notation of Dirac we write this expression as a “weak” relation: 


ġ ~ {9, H lan, (A.15) 


where Hr denotes the “total” Hamiltonian, Hr = H +um@m. The “~” means that one has 
to evaluate all Poisson brackets before setting the constraints to zero (Øm œ% 0). 
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Since the constraints Øm are functions of qn and pn as well, they must of course fulfill the 
same equation of motion (A.14) as the functions g. Hence, consistency demands: 


U bm x {om H Les + Um {Pm dm’ }pp : (A.16) 


This equation can be used to determine the um unless the second Poisson bracket vanishes 
({dm, m }pp=0). In that case one gets further constraints, so-called secondary constraints: 


x(q, HI © 0. (A.17) 


Of course the secondary constraints have to fulfill the equation of motion (A.14), too. This 
may lead to further secondary constraints and so on. 

Once all secondary constraints have been found, one can start to classify them. According 
to Dirac there are two “classes” of constraints: first class and second class (not to be 
confused with primary and secondary). A phase space function (or constraint) is called first 
class when its Poisson brackets with all other constraints is (weakly) zero. If at least one of 
these Poisson brackets is unequal zero, one speaks of a second class function/constraint. 

Let us go back to equation (A.16): As long as the second Poisson bracket does not vanish, 
one gets solutions for um: 


Um = Um(p, q) + Va Vam, (A.18) 


where Um are special solutions of the inhomogeneous equation and Vam are solutions of the 
homogeneous equation 


Vam{ bj, mtrs = 0. (A.19) 


The va are arbitrary parameters, which means that some kind of freedom is contained in 
the theory. Consider a dynamic variable g(t) with an initial value g(0) = go: After the 
infinitesimal time interval ôt one has 


g(5t) = go + got = go + {9, Hr}rndt , (A.20) 
where was used. Defining 
H < HL Ud and du S Vantin (A.21) 
(see {A.18)), one may write 
g(5t) = go + ôt ({9, H'a + Valg, Patre) - (A.22) 


Due to and the product rule for the Poisson bracket it is obvious that a is a first- 
class constrain{?} Furthermore, H” is also a first class function by construction (cf. (A.16)). 
As noticed earlier, va are arbitrary parameters, which means that g(dt) is ambiguous as 
well: Replacing vg with some v/, in leads to a different value for g(dt), the deviation 
being 


Ag(ôt) = Ôt (Va — vg ){g, Patre = €at9; Qa fre- (A.23) 


* £5, da Lp = Vam{j, om} PB + {¢;, Vam}pBdm = U. 
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If one interprets as a gauge transformation, then the first-class constraints a are 
obviously its generators. In doing two successive gauge transformations of g, one can easily 
show that the Poisson bracket 1 da: da: trn generates a gauge transformation as well. By 
applying the product rule one can furthermore show that the Poisson bracket of two first class 
constraints is a first-class constraint itself. Hence, {¢a, dar Ern must be a linear combination 
of the first-class constraints in the model under consideration. Therefore, we deduce that 
all primary and secondary first-class constraints generate gauge transformations} This fact 
should also be taken into account in the equations of motion. We therefore define the 
extended Hamiltonian 


Hg = Hr + Vida! . (A.24) 


The generators ġa’ are all those which are not already contained in Hr and are therefore 
first-class secondary constraints. The corresponding equations of motion are now given by 


ġ ~ {9, Hejes - (A.25) 


What about the second class constraints? In order to treat those we first consider the 
matrix Cap = {94, B Len where the 0.4 now denote all constraints, and for simplicity we 
assume the irreducible case, i.e. that all 64 ~ 0 are independent from each other. Obviously, 
det Cap & 0 if there is at least one first class constraint among the ġa. Redefining the 
constraints as 64 > a P Œn with an appropriate invertible matrix a F one can always find 
an equivalent description of the constraint surface in terms of constraints Ya % 0, Xa %0, 
whose Poisson bracket matrix reads weakly 


Ya Xa 
Yb U 0 
É E IE (A.26) 


where Cga is an antisymmetric matrix that is everywhere invertible on the constraint surface. 
In this representation, the constraints are completely split into first and second classes, and 
the number of second class constraints is obviously even. 

A possible way of treating the second class constraints was invented by Dirac in introduc- 
ing the so-called Dirac bracket 


(i: g9}os = LL dtes — 1, XatesC 3. g}rs, (A.27) 


where C? is the inverse of Cap. Since the extended Hamiltonian is first class, one can easily 
verify that H c still generates the correct equations of motion in terms of the Dirac bracket: 


ò ~ {9, HC Lan (A.28) 


The original Poisson bracket can be discarded after having served its purpose of distinguish- 
ing between first-class and second-class constraints and all the equations of the theory can 
now be formulated in terms of the Dirac bracket (see ref. [6] for a detailed proof). 


Initially, Qa = Vambm consisted only of primary constraints. 
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A.1.2. Field theoretic extension 


We are now interested in the field theoretic extension of the formalism developed above and 
illustrate this with an example: free Maxwell theory. The action is given by 


1 
Re / dt T dr Fa F, (A.29) 


with the electromagnetic field tensor F, defined in Eqn. (1.10). The fields A, (t, 71 corre- 
spond to the qn(t) in the previous section. The variable 7 can be interpreted as a “contin- 


uous” index. According to (A.6) with Ay = OA, = ou the conjugate momenta are given 
by 


j= (1 [a'r arra = F(a 
a) = FE s [E Feee) = FOE). (A.30) 


In analogy to {qn, Pn’ Ern = Onn’ We now have 
{A (2), r” (@)}pn = 646° (Z — 2’) . (A.31) 
Due to antisymmetry of the field tensor, equation yields the primary constraint 
T (Z) = 0. (A.32) 


A Legendre transformation, as defined in (A.5), gives us the Hamiltonian of Maxwell theory 
3 HÀ 1 rs 1 TU 
H = | x| nA, t EU Fps + 2P Fro l. (A.33) 


where Latin indices run from 1 to 3. The constraint (A.32), partial integration and the fact 
that Fro = =r" yields 


1 1 
H= [ee (Grek. + gut = oc . (A.34) 


All time derivatives have now been replaced by conjugate momenta enabling us to use the 


consistency condition to get 
Oxi? a {r?, Han = 0,0", 
which yields the secondary constraint 
Onn ~ 0. (A.35) 


A further consistency check shows that (A.32) and (A.35) are the only constraints, since 
{0,7", H ypg = 0. Furthermore, they are first-class because of 


{n° @) n° (2) jon = 0, {n° (2), O-n" (Z")}on = 0, 
{O,n" (Z), Opn” (Z')} pp = 0. 


Obviously, the Hamiltonian H is first-class as well and therefore can be used for H’ from 
(A .21)). The total Hamiltonian Hr hence becomes 


1 1 
Hr = f (Gren. + sel Pr - | Aoda" + [oE ds, (A.36) 
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where v(7) is arbitrary. Inserting Ag into the equation of motion (A.15), we see that 
v(#) = Ap(#). This means that the time derivative of Ag is ambiguous and that Ag as well 
as its conjugate momentum 7° = 0 are unphysical. Using the “extended” Hamiltonian Hg 
one may eliminate these unphysical quantities: 


Hpg = Hr + I (A.37) 


Choosing v(x) = 0 and u'(x) = u(x) — Ao one arrives at the new (simplified) Hamiltonian 


(cf. (A.36)) 


1 1 
H= i] (Gren. -+ Zara") dr + [voada (A.38) 


which still produces the correct equations of motion for all physically relevant variables. 


A.2. Natural Units 


In general it is always convenient to use a unit system where the formulae under consideration 
are rendered particularly simple, i.e. with only a minimal set of constants. In particle physics 
this is achieved by setting the velocity of light c and Planck’s constant h equal to 1 and 
dimensionless: 


c=h=1. (A.39) 


Quantities computed in these units can always be converted to SI units at a later point, if 
desired. However, let us take a closer look at the properties of these so-called natural units. 
From the speed of light being now dimensionless it follows immediately that 
[length] . 
1 = |d} = ——_— => [length] = [time]. A.40 

DE length] = [time (A.40) 
Time and length now have the same dimension. This also becomes clear when considering 
the relativistic formulation of space-time where the four-vector x” = (ct, 71 with c = 1. 

Similarly, one finds: 


1 


1 = [A] = [energy][time] = [time] = 


So time (and hence length) now has dimension of energy "T. In particle physics one obviously 
has to deal with the energy of particles, which is to small for the unit Joule to be practical. 
Therefore, one uses the unit of electron volts instead, where 


leV = 1,602177 - 10° 7, (A.42) 


which corresponds to the kinetic energy an unbounded electron gets when accelerated by 
an electrostatic potential of 1 Volt (Volt = Joule/Coulomb, cp. charge of an electron from 


Table[A.1). 


From E = mc? follows furthermore that mass now has dimension of energy as well, i.e. 


[mass] = [energy], (A.43) 
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and if we consider Boltzmann’s constant kg = 1, temperature also has the dimension of 


energy ({temperature] = 


=d T 
cuo’ 


[energy]). In the same way, uo = 1 and dimensionless leads to 
voltage and current having the same units (cp. Table A-1). From the relation €9 = 


it 


follows additionally, that the dielectric constant in vacuum en = 1 and dimensionless as well 


in natural units. 
Typical values are: 
10° 7eV: 
some eV: 
5,11-10°eV = 511lkeV: 
some MeV (i.e. 10°eV): 
938MeV = 0,938GeV: 


speed of light: 


thermal energy of a particle at room temperature 


chemical reactions 

rest mass of an electron 
nuclear reactions 

rest mass of a proton 


c = 299792458” 


Planck’s constant: 


unit charge: 
gravitational constant: 


h = 6, 62607 - 10-34.Js 


h = Æ =1,05457- 10734 Js 


27 
e = 1,602177-10-"C 


_ —11m3 
G = 6.673 - 10 pore 


(electron: 


Boltzmann’s constant: | kp = 1, 3807 - 10-845 
uo = 4r - 107 


magnetic constant: 


Table A.1.: some constants of nature 


A.3. Dimensional regularization 


We follow the derivation given in reference [3] on pages 382-385: Working in n-dimensional 
Minkowski space we start with integrals of the type 


1a) = [oa (A.44) 
nid ~~ (k2 + 2kq — L?)*’ è 
where k = (ko,r,ġ,01, 02, ...,On—-3) in polar coordinates and the volume element is therefore 
given by 
n—3 l 
d"k = dkor” *drdé ] | sin’ 0:00. (A.45) 
i=1 


Shifting variables k/, = ky + q and using 


n—1 


20 n—-3 T 
f do J] I sin’ 0,40, = i, (A.46) 
S i=1 9 EES) 


2 


we arrive at 


n—1 co co 
2m 2 r?— 2dr 
In = - pa Ts A.A7 
OTe LO Ge = 


where the primes have been dropped. Since the integrand depends only quadratically on 


ko we can replace T dko > 2 i dko. To evaluate the remaining integrals we use the Euler 
—oo 
beta function 


T R(x) > 0 
B = = 2 |] att se AA 
ay) pee a2 l ape I. AER (A48) 
0 
Considering 
1+ 1+ 8 s 
= — = Qa — —— t= — A.4 
pa yaa i, (A.49) 
leads to 
F sÊ L os Dla— = 
, (A.50) 


|e + M?)e = 2 LM UBS (a) 


and if we identify 8 = 0 and M? = —r? — (q? + L?) we can use this formula to perform the 

integration over kg in S 

[ad HHT 
(P+ TAT 


Tn(q) (A.51) 


0 


By identifying a’ = a — 3, B = n — 2 and M? = q? + L? we can once again use formula 
{A.50) to perform the remaining integral: 


= = Cm De) ge =a (A.52) 
Therefore the result is 
dk 2. 2 2 dud 
rlo = l —— = (- + (- 94° — L oe ; 


Differentiating both sides with respect to q, and redefining a as well as using the property 
of the Gamma-function zT (x) = T (1 + 7x) leads to 


h(a) = | Perr = h), (A.53b) 


and further differentiation with respect to qy yields 


kťk” g” (—¢? = I?) 
(= | kee = | d + |] Lg). A. 


A useful trick in order to bring integrals appearing in typical Feynman graphs into the 
form I,(0) or If"(0) as defined above, is given by Feynman’s formula 


1 
1 dz 
ab J (az + b(1— 2)? Gia 


The proof of this integral formula is straightforward: simply substitute z’ = (a — b)z — b. 


Al 
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